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Abstract. We resolve the 'baxterization' problem with the help of the au- 
tomorphism group of the Yang-Baxter (resp. star-triangle, tetrahedron, . . . ) 
equations. This infinite group of symmetries is realized as a non- linear (bira- 
tional) Coxeter group acting on matrices, and exists as such, beyond the narrow 
context of strict integrability. It yields among other things an unexpected elliptic 
paramctrization of the non-integrable sixteen- vertex model. It provides us with 
J> ■ a class of discrete dynamical systems, and we address some related problems, 

\ such as characterizing the complexity of iterations. 
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1 Introduction 



The results presented here originate in a long-standing collaboration of M. Bel- 
Ion, J.M. Maillard and the author [jl], |, g, |, g, g, |], augmented by further 
elaborations (as || ||). 

The Yang-Baxter equations and their variations (star-triangle equations, 
tetrahedron equations,...) are nowadays considered as a characteristic feature 
of integrability of mechanical systems (cla ssical hamiltonian systems, as well as 
quantum and statistical systems) JToJ y, || [l| 0, H jl|, |o) . . . . 

This field has expanded very much for the last twenty years. We shall not dwell 
here on its innumerable developments, and relegate general considerations to 
the end. 

We shall start from the Yang-Baxter equations and view them as an algebraic 
system of equations on matrix entries. 

We first recall (section ||) what the two most common forms of the equations 
for vertex models are. 

In section || we introduce certain groups generated by involutions (Coxeter 
groups) [ pT| , together with some realizations in terms of birational transforma- 
tions of projective spaces. These realizations are generically denoted T in the 
sequel. 

In section || we build up, from the groups T, an infinite group of automor- 
phisms (denoted A) of the Yang-Baxter equations. One should keep in mind 
that the equations we analyze form an over deter mined system, making the ex- 
istence of a large group of symmetry quite remarkable. 

We next show, for the paradigmatic example of the Baxter symmetric eight- 
vertex model and with both a picture and algebraic results, how our automor- 
phism group reconstructs the well known elliptic curves of solutions. This is a 
first solution to the 'baxterization' problem J22|. 

In section || we describe more groups of matrix transformations. These 
are very similar to the ones of section |[ but are adapted to the star-triangle 
relations, tetrahedron equations, and higher dimensional generalizations. 

In section we complete the solution to the baxterization problem, empha- 
sizing the notions of T-covariant versus T-invariant quantities. 

We next produce (section an algorithm to calculate the covariant polyno- 
mials as well as the algebraic invariants. 

One of the features of T is that it is defined outside of the space of so- 
lutions of the Yang-Baxter equations, and we may look at its action beyond 
integrability. It becomes an ordinary discrete dynamical system on the space 
of parameters, possibly non-conservative. In section ^|, we describe this action 
for the general 16-vertex model, i.e. the most general two-state vertex model 
on a two-dimensional square lattice, which is known not to be integrable. The 
outcome is quite surprising: the orbits of T stay within 1-dimcnsional subvari- 
eties of the 15-dimensional space of parameters. This reveals the actuality of an 
amazingly large number of algebraically independent invariants of T, and yields 
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an elliptic parametrization of the non-integrable model, without reference to the 
Yang-Baxter equations. 

In section |ll| we turn to the notion of complexity of the iterations and show 
how it fits with the existence of algebraic invariants. 

In the last section (Q) , we exemplify some of our results on a definite system, 
related to a three dimensional vertex model. 



2 The Yang-Baxter equations (vertex models) 

Many presentations of the Yang-Baxter equations appear in the literature, di- 
viding into two classes: 

The first class contains a parameter, also called spectral parameter |Q, and 
reads, for vertex models: 

J2ai,a 2 ,a 3 R ala 2 (^li ^) R jia 3 ( A l > X ^) R hjT ' CO 

= E/9i,ft,,ft ^| 3 (A2,A 3 ) J R^(Ai,A3)-Rf 1 1 / 2 2 (Ai ) A2). 

In this system of equations, R is a matrix of size q 2 x q 2 , whose entries are 
functions of the parameters Ai, usually through the difference Xi — Xj (at least 
when A labels a point on an elliptic curve). 

The parameter first appeared as parametrizing the solution. It was later 
understood as a "spectral parameter" in the context of the quantum inverse 
scattering method. The origin of the parameter is one of the issues we wish to 
address. 

A given family of models is obtained when one specifies the size q and the 
relations between the entries of R. Typical relations are equalities between 
different entries or vanishing of some others (see for example section ||) . 

Since the entries of R are allowed to be different for the three copies of R 
entering (0), one may rewrite (|l]) with three explicitly different matrices A,B, 
and C of the same size q 2 x q 2 . 

EAiii 2 R"l«3 f t>2t«3 _ \ " /^i2«3 T3il03 A0102 (n\ 
/i aia2- D iia 3 j2j3 ~ ^0203 ° Pi 33 A hh \ Z > 

"1,02,03 Pl,02,/33 

or shortly 

A12B13C23 = C23B13A12 (3) 

with now usual notations. The latter form has the interest of not referring to 
any explicit parametrization. 

The second class contains no parameter, and is sometimes called "constant" 
Yang-Baxter equation. It reads: 

i?i2-Ri3^?23 = R23R13R12 (4) 
that is to say like (|) but with A = B = C = R. 
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To go from equation ([!]) to equation (^) requires essentially an adequate 
choice of Ai, A2, A3, namely Ai = A2 = A3. The reverse move, i.e. to recover the 
parameter dependence from a constant solution was given the name of 'baxter- 
ization' p2| . 

The use of the spectral parameter is clear in the context of the Bethe ansatz, 
since it serves building up the generating functional for commuting conserved 
quantities, but we shall not be concerned with this aspect. Its origin will be 
explained in the next sections. 

3 Some operations on matrices: the groups T 

We describe here some elementary operations on matrices of various sizes. The 
matrices wc consider arc defined up to an overall multiplicative factor. The 
space of parameters is thus some projective space P„ with the entries of the 
matrix as homogeneous coordinates. 

Let us start with the matrix R of a two-dimensional g-state vertex model on 
a square lattice. R has the structure of a tensor product and the indices of R 
are pairs of indices: 

R u i,j,k,l = l..q (5) 
We may define the inverse I up to a factor (well defined in the projective space) 

i,j,k,l = l..q (6) 

a/3 

with fi an arbitrary multiplicative factor, the transposition t: 

{tR)l = R?i i,j,k,l = l..q (7) 
as well as two partial transpositions ti (index I for left) and t r respectively by: 
(tlR)%=Ru (t r R)%=R% i,j,k,l = l..q (8) 

Of course 

t = tit r = t r t h I 2 = t 2 = tf = t 2 r = 1, and It = tl (9) 

However 

U I ^ I U and t r I ^ I t r (10) 

The two partial transpositions do not commute with the inversion while their 
product t does. The transformation til and t r I are of infinite order. Notice 
of course that they may have finite orbits when acting on certain non-generic 
matrices. 

The group T generated by i; , t r and / essentially consists of the iterates of 
til and its inverse Iti (or t r I and It r ) up to multiplication by an element of 
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its center. It is generated by involutions realized as rational (possibly linear) 
transformations of some projective space P„. 

We shall describe in section [| more instances of groups T acting as rational 
transformations of matrices, and generated by involutions. 



4 The automorphisms of the Yang-Baxter equa- 
tions 

Suppose we have a solution (A, B,C) of (||): 

Al 2 Bi 3 C 23 = C 23 Bi 3 Ai 2 (H) 

We may take the partial transpose t\ of ([n]) along space 1 and get: 

(tiBia) {hA 12 ) C 23 = C 23 (hA 12 ) (t x B 13 ) (12) 



Taking the partial transpose t 2 of (12) yields 

(tiSia) (t 2 C 23 ) (ht 2 A 12 ) = {ht 2 A 12 ) {t 2 C 23 ) (hB 13 ) (13) 
that is to say: 

{tiB) 13 {t t C) 23 {tA) 12 = {tA) 12 (t t C) 23 (t t B) 13 (14) 

Multiplying both sides of ([[ij) to the left and the right by the inverse It A of tA 
gives: 

(tIA) 12 (t,B)i3 (UC) 23 = {UC) 23 {UB) l3 (tIA) 12 (15) 
which shows that the transformation 

K A : {A, B, C) — > (tJA, 4,5, t,C) (16) 

takes a solution of (||) into a solution of @. 

We could similarly construct Kb (resp. Kc) which inverts and transposes 
B (resp. C) and acts by partial transpositions on A and C (resp. A and £>). 

The automorphisms Ka,Kb,Kc are involutive and generate an infinite group 
of automorphisms A of the Yang-Baxter equations. It is possible to show that 
A is isomorphic to the Weyl group of an affine Lie algebra of type A% || 0. 

The group T appears when we look at the action of A on one of the indi- 
vidual copies of R entering equation (Q). However, the action of T is defined 
without reference to (Q), and we will now concentrate on it, and somehow forget 
about (lil). 
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5 Baxterization of the Baxter model 



The exemplary two-dimensional integrable vertex model is the Baxter model [ p3| , 
p4| , tt2j, with matrix of Boltzmann weights: 

(a d\ 
\d a/ 

where one sees the vanishing conditions and equality relations between the en- 
tries. It is straightforward to see that the form (JIJ) is stable by all the generators 
of r. We say that the pattern (p7|), i.e. the collection of vanishing conditions 
and equalities between entries is admissible (see [|[ ||). 
The transformations t\ and t r read: 



b, c — ► d, d —>■ c, 



and / reads 



a* 



d? 1 b 2 -c 2 ' 6 2 -c 2 ' a 2 -d 2 ' 



One may look at the actiofr-ofT on R. Take one starting point in the 
space P3 of parameters, iterate on it -and just look at the orbit. One gets 
the following picture, showing a three-dimensional perspective of the orbit, in 
inhomogeneous coordinates u — b/a, v — \c/a, w = d/d, with starting point (*): 



Figure 1: 'Baxterization' of the Baxter model: a perspective view 

One sees from this figure that the orbit lies inside a curve. The action of Y 
must have at least two algebraically independent invariants. Indeed there are 
two such invariants: 

o 2 + b 2 — c 2 — d 2 ab — cd 

' ab + cd 2 ab + cd 
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The curves Ai = constant, A2 = constant, are precisely the elliptic curves 
which appear in the solution of (jl]) , but we have obtained them with the form ( fi"| ) 
as only input, and in particular without having to solve ([1J). The resolution 
of (|I|) becomes straightforward once the adequate parametrization of the above 
curves is used. The orbit we get, although discrete, tends to be dense inside one 
dimensional curves. The phenomenon is similar to the iteration of a rotation 
(obtained as the composition of two axial symmetries) with irrational rotation 
number. 

Warning: It is in no way the existence of these curves which ensures the 
compatibility of the Yang-Baxter equations! 



6 Other transformation groups 

There are many other equations like (|l]). They also have automorphism groups 
generated by (bi)rational transformations. This means we may construct many 
different realizations T. Here are a few more. 



6.1 Prom spin model with interactions along the edges 

As was shown in ||, the star-triangle equation also has an infinite group of 
automorphisms. The matrices under consideration are ordinary q x q matrices 
m, of entries m^. The generators of T are the matrix inverse / as previously 

*^2(Im) ia m aj = fi Sij (19) 

a 

with u an arbitrary multiplicative constant, and the element-by-element inverse 
(Hadamard inverse): 

J : mij > l/rriij (20) 

/ and J are two non-commuting rational involutions. The product ip — I J is a 
birational transformation of infinite order. 



6.2 Prom vertex models in three (and more) dimensions 

The matrices under consideration in this case again have multi-indices. The 
matrices are of size q 3 x q 3 , with entries of the form: 

RlL (21) 

By examining the symmetries [EL || of the tetrahedron equations [^5|, ^6|, |27| : 

-^123^145^246-^356 = -^356-^246-^145-^123- (22) 



G 



with the usual notations, one is lead to the group T generated by the inverse I 
and the three partial transpositions t g , t m and ts, with 

and similar definitions for t m and i^. 

One could as well consider multi-index matrices of size q d x q d written in 
the form M' 1 ' 2 '"' 1 '. Section El deals with d = 2. The tetrahedron equation 
corresponds to d = 3, and so on. There exist d different partial transpositions 
ti, t%, . . . , td with obvious definitions. Combined with the inverse /, they gen- 
erate an infinite group T. We will return to the question of the size of these 
groups in sections || and 11 . 

6.3 More groups 

One may further enrich the representations by imposing constraints on the en- 
tries of the matrices, provided the transformations are compatible with these 
constraints (see |ij for the notion of admissible patterns). This yields realiza- 
tions on projective spaces of lower dimensions. 

If the dimension n of the projective space is 1, the only birational transfor- 
mations are homographies. If n = 2, the group, called Cremona group [p8| , is 
already much larger. Its elements may however be written as products of linear 
transformations and the basic blow-up J. The interesting subgroup made out 
of polynomial and polynomially invcrtible transformations of the plane has a 
more elementary structure [^9[ [}(J. When n is larger than 2, the structure of 
the group of birational transformations is much more complex. 

Let us stress that, at the level of the realization, there may exist additional 
relations between the generators, possibly making it finite. 



7 Baxterization 

It appears clearly in sections [| and || that, for the Baxter model, the existence of 
the parameter comes from the existence of sufficiently many algebraic invariants 
of r to confine its orbits, and therefore the ones of A, to curves. The general 
situation is more subtle, but it always amounts to the following: find the T- 
covariant varieties passing through a constant solution, and/or look for solutions 
only on covariant varieties. 

Different situations may appear, and we illustrate them by specific examples: 
- the simplest one was encountered in section ||, recalling that the only 
input information was the form of the matrix (|l7|). The parameter is just the 
uniformizing parameter for the elliptic curve which pops out of the action of T 
on any of its generic points. In this case, the action of T has INVARIANTS to 
imprison the orbit, and the entire space of parameters is foliated by invariant 
subvarieties (curves for the Baxter model). Each generic spectral curve has an 
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infinite group of automorphisms, and they consequently all have genus or 1. 
The solution of (||) is obtained in the following way: choose any point A in 
the space of parameters. Through A passes one curve of the foliation. Choose 
B anywhere on this curve. The last matrix C is completely determined then. 
What is special about the Baxter model is that there is no restriction on the 
first choice (A), within the space P3 of parameters. This is far from being the 
general occurrence. 

- a more subtle situation is encountered for the chiral Potts model (as well 
as for the free fermion condition pi], [32] |33|). Integrability appears on a alge- 
braic variety of which all points have a periodic orbit under T. This variety is 
automatically globally invariant by T, but there may very well be NO INVARI- 
ANT of the action of T. Besides, the resolution of the Yang-Baxter equations 
may lead to consider further algebraic conditions. Solutions are to be found 
on a possibly isolated subvariety of the space of parameters. In this situation, 
the spectral curve may have genus higher than or 1, since it has only a finite 
number of automorphisms |34|, |35|, p_9j. 

- a similar situation is encountered in the case of the Jaeger-Higman-Sims 
model |37|, 38 1 . The parameter space of the model is P2. The matrix of Boltz- 



mann weights belongs to a 3-dimensional abelian algebra of matrices of size 
100 x 100, with three generators \,A,J — 1 — A where 1 is the unit matrix, 
A the adjacency matrix of the Higman-Sims graph, and J the matrix with all 
entries equal to 1. Two products exist in this algebra, respectively the usual 
matrix product and the element by element product. The symmetry group is 
generated by the inverses / and J, corresponding to these two products, as in 
section (3.1), equations (pi, I2T 



C (linear map of P2) 



The two inverses are related by a collineation 

In terms of homogeneous coordinates [xq, x\, X2]' 

J : [xo,xi,xz] — ► [0:1X2, a; a;2,xiX2] 
"1 22 77 
C= 1-8 7 
1 2 -3 

We can find three T-invariant subvarieties of P2: one line D, and two hyperbolae 
HI, H2: 



D : 
HI 
H2 



X 2 - Xi 







X,2 + 3 Xi X2 — 3 Xo X2 — Xo Xi =0 
Xt X% + 2 Xi 2 — 2 Xq X2 — Xg Xi =0 



Hyperbola HI is the one of reference |39 . The two hyperbolae contain an 



infinite number of singular points, and we know this prevents the existence of 
an invariant but not of covariant polynomials, and a fortiori not of covariant 
ideals. 
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8 Covariants and invariants 



It is fortunately possible to go further in the analysis of the invariant subvari- 
eties |J. 

We start from a group G generated by by v involutions I\, I2, ■ ■ ■ , ifc, (k = 
1 . . . 1/), verifying no relations other than the involution property. The group G 
is infinite and there are two essentially different situations. 

If v = 2, the group is the infinite dihedral group Z2 x Z, and all elements 
may uniquely be written If (ii/2) 9 , with a — 0, 1 and q 6 Z. The number of 
elements of given length I is 2. 

If v > 3, the number of elements of length / grows exponentially with I, and 
the group is in a sense bigger (still countable). 



As an example, for the groups described in section 6.2, the number v of 
generators depends on the dimension d of the lattice: it is just 2 d ~ 1 so that if 
d = 2, G is generated by two involutions and if d > 3, G is generated by more 
than three involutions. 

We then construct various realizations T of G by explicit transformations 
of some projective space. They are obtained by specifying the realization of 
the generators. We use the same notation Ik for the generators of G and their 
representatives in T. The realizations r of G may be written as polynomial 
transformations in terms of the homogeneous coordinates. 

Each involution defines a characteristic polynomial (f>k of degree d\ — 1 in 
the following manner. The Ik being involutions, if appears as the multiplication 
by a degree d\ — \ polynomial <fik(xo, ■ ■ ■ ,x n ). 

A F-covariant polynomial P verifies: 



P(7(as)) = a(7,a;)P(a!) V 7 e T,Vx e P„ (24) 

The coefficient 0(7, x) has to fulfill the cocycle condition: 

0(7172)2?) = a(7i,72 2)0(73, x) (25) 

Indeed relation (|^) demands 

P{l\l2x) = 0(7172, x)P(ar) 

= a(7i,72 x)P(j 2 x) 

= a(7i)7a x)a(j2,x)P(x) 

The cocycle a will be completely determined by the values of a(Ik, x), k = 
These values may be found easily: when applied to 71 = 72 = ^fe, 
condition ( |25| ) shows that a(Ik,x) has to divide a suitable power of 4>k- 

Finding an invariant A = P/Q implies finding two polynomials P and Q of 
the same degree, which transform the same way (are covariant with the same 
cocycle) under all the generators, i.e: 

P(Ik(x))=a{I k ,x)-P(x) and Q(h(x)) = o(J fc , x) ■ Q(x) (26) 
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Once the cocycle a is chosen, solving (|2(j) becomes a handable linear problem, 
of which the compatibility can be further studied Q . We have proved that the 
proliferation of singularities impeaches the existence of any invariants, but the 
converse is not true. One of the outcomes is that, in this general setting, the 
existence of invariants is exceptional. 

Finding all invariant subvarieties is more tricky. What we have obtained 
so far is invariant subvarieties of codimension 1, or subvarieties determined by 
invariants of T. Smaller subvarieties have more than one equation, and (|24|) has 
to be replaced by a matrix relation: 

11(7(3;)) = Aft, x)n(x) v 7 e r, Vx e p„ (27) 

where II is a vector constructed from the (non canonical) list of equations, and 
A a matrix. In other words we demand that the ideal defining the subvariety be 
invariant by T, keeping in mind that the subvariety we look for does not have 
to be a complete intersection. 

Let us describe an example, coming from the hard hexagon model jlO], ftl[ . 
Consider the two involutions of P4 given in terms of homogeneous coordinates 
[xo,xi,xg,xs,X4]: 

h - [x eXl x 4 ,x 4 (x x 5 ~x 2 2 ) 1 ~x 2 x 1 x 4 ,x 1 {x x s -x t 2 ),xox lX4 \ 
* 2 — [x 4 x 2 x 5 ,- Xl x 2 x 5 ,x 5 (x x 4 -x 1 2 ),xox s x 5 ,x e (XOXJ-X! 2 )} 

The two inverses are related by a linear transformation, i.e ii = r i\r with 

r : [xo,x 1 ,x 2 ,x 4 ,x 5 } — > [x , x 2 , xi , x 5 , x 4 ] (28) 
and 4>\ = i\ = xj 2 x 4 2 (xo xg — x 2 2 ) 2 ■ If 

III = Xq 2 X2 2 X4 + Xq 2 Xi 2 Xs — Xo 3 X4 £5 — X4 2 X§ X-} 2 — X4 X^ 2 Xl 2 + Xt^X^ 2 — Xl 2 X2 2 Xo 

II2 = xo 2 — X4 x 5 — A xi x 2 
with A a free parameter, then we do have relation ( p^ ) with a matrix 

At = 



X! 3 X 4 3 (x x 5 - x 2 2 ) 

-x 4 ~x 4 2 (x x 5 - x 2 2 ) 



and 



Ax(ix{x))Ai{x) 



4>\(x) 
4> 2 (x) 



Both III and II2 are unchanged by r, III is covariant, and consequently the 
variety III = is invariant by T. In contrast, II2 is not covariant and only the 
intersection of 1I2 = with III = is invariant by T. 
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9 The sixteen-vertex model 



We abandon here the strict context of integrability, and turn to problems of 
discrete dynamical systems, defined in the space of parameters of the models. 

The field of discrete dynamical system is a well developed one. The notion of 
integrable map goes back, in the context of hamiltonian dynamics, to Poincare, 
who founded the subject H|, see also |44[ |||]). The last 30 years have 
seen the subject expand a lot, specifically with the advent of computer calcula- 
tions. In order to gain simplicity, one was lead long ago to studying iterations 
of polynomial and polynomially invertiblc transformations of the plane fltl ), 
eventually renouncing to hamiltonian structures. Analytic maps have also been 
fruitfully analyzed, especially in one complex dimension |47], Q . More recently, 
some remarkable multi-parameter families of maps, eventually integrable, have 
been constructed from soliton equations |5^, |5l], |5^, |5^, Q . 

We want to analyze the behaviour of iterations of a typical infinite order 
element of our groups T such as the product of two generating involutions. The 
first example we will examine is of great interest for statistical mechanics on 
the lattice. It comes from the general 2-state vertex model on the square lattice 
in dimension 2, i.e. the sixteen-vertex model. Take R to be the general 2x2 
matrix. This matrix up to a multiplicative constant represents an element of the 
15-dimensional projective space P15. A typical orbit of the adequate T, when 
projected to a 2-plane of coordinates, looks like: 




Figure 2: One generic orbit ofT for the 16-vertex model 
This figure shows how efficient the graphical method can be. The situation is 
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very favorable since the orbit is confined to a low dimensional variety (a curve). 
What can be proved B, following the ideas of section (^), is the existence of 
a collection of invariants, of algebraic rank 14. More precisely, there are 17 
invariants with 3 relations. As in Figure 1, the discrete orbits are generically 
dense in a curve, and this helps making the 'graphical detector of invariants' so 
useful. 



10 More pictures 



One should not believe that T-orbits are always curves. Let us introduce a class 
of transformations in P2 which contains in particular the transformations (pi]). 
Consider T generated by I and J, given in homogeneous coordinates [xq, X\, x%\ 
by: 



J : [x ,xi,x 2 ] — > [xix 2 ,x X2,xix 2 } 



(29) 
(30) 



with C a projective linear transformation of P2. This class of realizations, 
parametrized by the collineation matrix C already contains all kinds of be- 
haviours. Moreover, many nice examples obtained from matrix inversions do 
implement relation (|3C 



10.1 Non-symmetric Z 7 

Suppose m is the following cyclic 7x7 matrix: 



The matrix inverse verifies (0 with| 



X 


y 


y 


z 


y 


z 


z \ 


z 


X 


y 


y 


z 


y 


z 


z 


z 


X 


y 


y 


z 


y 


y 


z 


z 


X 


y 


y 


z 


z 


y 


z 


z 


X 


y 


y 


IS 


z 


y 


z 


z 


X 


X / 


y 


y 


z 


y 


z 


z 



Cz 7 



2 6 6 
2 -l-n/7 -l + n/7 
2 -i + i-v/7 -l-n/7 



(31) 



(32) 



The iteration of ip = IJ, with / the matrix inversion and J the element by 
element inversion, yields the following picture, in the inhomogeneous coordinates 
u = y/x, v = z/x, if we represent a number of orbits at the same time: 

2 This form is a particular value of an eigenmatrix of conference digraph [ jssfr 



12 




Figure 3: A collection of chaotic orbits 

Figure 3 shows there should not be any invariant. Indeed, as proved in J9|, 
there cannot be any because T has an infinite number of singular points. The 
first singular points are visible on the figure. Their location is very easy to 
compute. 



10.2 A finite diagram model 



Take 



C 



2 
1 1 

-1 1 



(33) 



This has been introduced in || and the denomination 'finite diagram model' 
refers to the finite number of singular points of T. Although such a transforma- 
tion was not directly obtained from a matrix inverse, it (or similar ones) appears 
in properly defined reductions of the inversion of matrices of larger sizes p6| . 
We have found no invariant for this model. The graphical method is quite inef- 
fectual at generic points. It proves nevertheless useful in the vicinity of a fixed 
point. 
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-q. 




Figure 4 



Figure 4 represents a collection of orbits not far away from a fixed point 
of (/J) 2 . The fixed point [a 2 , a, — 1], a = (1 + \/5)/2, {x\/xq ~ .618, x 2 /x ~ 
— .382) is an elliptic fixed point at which the eigenvalues of the linear differential 
are the roots of s 2 — 17/16 s + 1. We see that the orbits explode when one gets 
away from the fixed point. This is exactly the situation encountered in J57J, 
see M. 



11 The complexity of iterations 

More insight into the behaviour of the iterations can be gained by looking at 
their complexity, in a way we will define. 

The first grasp at the complexity of T is through its size (section |^) , but we 
shall not comment on this here. 

We may satisfy ourselves with examining the simple case v — 2, where the 
group is essentially two copies of Z, and describe in the next section a specific 
example (section [12]) with v = 3. 

The simplest situation is the one where T acts by transformations of the 
2-dimensional projective space P 2 - We have given in section [To] a class of such 
a realization, and the study of 'the dynamics of singularities' as explored in |^| 
proves extremely useful. 

It is nevertheless something else we have in mind here: suppose ip = IJ is of 
degree w. Its iterates ip k do not have to be of degree w k . Indeed there is a simple 
mechanism for the lowering of the degree d(k) of ip k , for one has to factorize out 
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common factors from the expressions of the homogeneous coordinates of ip^ . 
We may analyze the sequence d(k), as a function of the order of the iteration 
k. This measure of the complexity of the iterations actually coincides with the 
one introduced by Arnold (5^] ; see |6(J for the particular case of bi-polynomial 
transformations of the plane. It brings a variety of behaviours, lying between 
the generic exponential growth and periodicity, with the particular instance of 
polynomial (or polynomially bounded) growth. In a given system of coordinates 
the notion of degree is well defined. Of course, changes of coordinates may 
change the degree, but allowed birational changes will preserve the nature of 
the growth. 

We have calculated the degree of the successive iterates of <p in a number of 
cases, the simplest ones being given by the class ( p9| , ^0| ). 
If C = Cz 7 given by (B2|), the first terms are: 



1- 



4, 12, 



33, 



232, 



609, 



(34) 



This may be seen as the first terms of the sequence: 

k 



i=0 



with fi the Fibonacci sequence: 
1, 2, 3, 5, 



13, 21, 34, 



This behaviour is found for any generic elements of the four parameter fam- 
ily |3l) of collineations, 



1 p q — p — 1 
Is -s - 1 
1 r —r — 1 



(35) 



Notice that the Jaeger-Higman-Sims model belongs to this family (compare 
collineations ( p4] ) and (pq)), and has the same sequence (p4|). 

For the 'finite diagram model' with C given by ( |33| ) the degree reaches its 
maximal value: 



d{k) 



For 



C{w) 



1 w-1 
1 -1 
1 



w 


-1 



(36) 



r is known to have an invariant, and the orbits lie within elliptic curves of 
equation 



P^(x , xi,x 2 )Q w (x , xi,x 2 ) - X(xx + x 2 ) 4 (x - x 2 ) 2 (x - Xx) = 
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where 

Pw = (1 - w)(x\ - XQXl) + (w - 3)x 2 (xo - Xl), 
Qw = (1 - w 2 )(x^ — XoX 2 ) + (w 2 — 4w — l)a; 2 (a;o — x 2 ) + 2(u> — l) 2 2i2;2(£o — ^1) 
and the sequence of degrees starts with: 

1, 4, 10, 20, 34, 53, 75, 102, 132, 167, 206, 249, 295, 347, 402, 461, . . . 

that is to say the second derivative d(k + 1) — d(k — 1) — 2d(k) is periodic of 
period 12. Notice that ( |36|) has the form (|35|). It is just not generic. 

As a last example, in the Baxter model, and also for the 16-vertex model, 
the iteration of <p = t g I yields the sequence of degrees: 

1, 3, 9, 19, 33, 51, 73, 99, 129, 163, 201, ... (37) 

from which one infers 

d(k) = 2k 2 + 1 

We have examined many more examples, and they all lead to the conjec- 
ture chaotic behaviour implies exponential growth, as regular behaviour im- 
plies polynomial bounds. 

The proof of this conjecture should follow from considerations on addition 
on elliptic curves f^f at least for realizations in P2 , the main idea being that <p 
is a constant shift on some curve. 

Another approach goes through the study of the factorization properties: the 
reduction of the degree d(k) from exponential to polynomially bounded comes 
from a sufficiently regular factorization process. An analysis of this process was 
successfully undertaken in J63| , |56"| for various realizations T and confirms up to 
now our conjecture. We shall give one more example in the next section. 



12 An example 



In J7j we introduced a restriction of the general 2-state model on the cubic lattice 
in three dimensions, by imposing the following relations on the entries of the 
i?-matrix: 

RfXl - if ' 1 '•.".'../:./•_'./.'. 1 (39) 

R is the direct product of two times the same 4x4 submatrix ||. It is further 
possible to impose that this 4x4 matrix is symmetric, since the product t = 
tit m t r acts as its transposition. 

J1J2J3 i Hll 2 »3 
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We use the following notations for the 10 homogeneous entries of this 4x4 
submatrix: 



( a 

di 
d 2 
V d 3 



dx 
bi 

(-2 



d 2 

C3 
h 



d 3 \ 

C 2 
Cl 
h J 



(41) 



The generators of T are the inversion of R, which acts as a matrix inverse on 
(fll]), and the permutations of entries ti, t m , and t r being respectively: 



ti : 


C3 


^ <^3, 


Cl 


<-> <i 2 




Cl 


<-> di, 




<-> d 3 




C2 


<-> d 2 , 


Cl 


d\ 



A numerical iteration gives the following picture of an orbit, projected on a 
2-plane of coordinates. 



























V- ■ 


^ ^^^^ 


w 



Figure 5: A 'two-dimensional' orbit 



This is another example of the efficiency of the graphical detector: we see 
that the orbit is confined to a 2-dimensional surface. This may be easily proven. 
If one defines 

P3 = ab 3 + bib 2 - c\ - d\, q 3 = c\d\ - c 2 d 2l (42) 



17 



together with the polynomial p 2 and p\ obtained by permutation of the indices 
1,2, and 3, and (with a corrected misprint from 0): 

r 3 = ab 3 - &1&2 + c 2 - d§, 
S3 = (a + b 3 )c 3 - did 2 - cic 2 , 
f 3 = (6 X + b 2 )d 3 - d x c 2 - C!d 2 - 

then 

• The polynomials pi, qi, i = 1, 2,3 form a five dimensional space, and any 
ratio of these polynomials is invariant by the whole group T. This means 
that the orbit of T is confined to a 5-dimensional subspace of the original 
9-dimensional parameter space. 

• If we consider the subgroup generated by / and t r , then r 3 ,S3,t3 are 
three more covariant polynomials, and they furnish three more invariants, 
proving that the orbit is confined to 2 dimensions. 

As far as the complexity is concerned, the non-chaotic image fits with the 
sequence of degrees, which is the same as (|37]). The analysis of the degrees 
can be performed easily: Let M be a generic matrix of the form (pLT|), and 
let Mo = t r M be the starting point of the iteration. We iterate t r I defined 
polynomially in terms of the entries, and keep track of the common factors 
(factorized out from the result of the action of K). Let Mk+i = t r I(Mk), ft 
the extracted factor, and = det(Mk). One gets do = det(t r M), fa = 1, 
d\ = det(t r It r M), f 2 = 1. Then a factorization of the determinants appears 
regularly, defining in the course an additional sequence Sk such that: 

A 2 = Af,S 2 A 3 = A? S 3 A 4 = 5|J 4 A 5 = Sj6 5 ... A n = Sl_ 2 S, 
/ 3 = A§ U = Aj h = AH A% = 81 ... f n = Sl_ 3 

Denoting by Uk,Vk,gk, Xk the degrees of respectively A^, Sk, fk and of the matrix 
elements of in terms of the entries of Mo, this implies 

u k = 3 Vk-i + v k 
9k = 2 v k - 3 
u k = 4n 
x k = 3 x k ~i - gk 

resolved by 

Xfc+3 - 3 x k + 2 + 3 x k +i - x k = and d(k) = 2k 2 + 1 

This is an instance of the general result obtained in j56) (formula (3.19), for 
q = 4), giving a good flavour of the method. 
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13 Conclusion and perspectives 



Inversion relations have proved to be a powerful tool in the study of statistical 
mechanical models, leading in particular to exact relations on the partition 
functions, even for non-integrable or higher dimensional models 65, jocf . 

The description often appeals to a diagrammatic representation 67 ] , which 
we have not used here. 

The Yang-Baxter equations in their various forms have brought an enormous 
amount of exact results in the field. They have unified branches of physics and 
mathematics such as ice models [|8[ ^] and solitonic wave equations ]70| fn}, 
gave new insights into knot polynomials fl% , |73| , f74[ f75f , and produced successful 
offsprings like quantum groups j[6], [78| . 

What H contain is a concrete link between discrete dynamical systems 
and quantum integrability, producing a large number of interesting dynamical 
systems. We have been able to use a little part of the paraphernalia available 
from the field of dynamical systems -the simplest being numerical calculations-. 
One outcome is that the rational nature of the transformations captures a good 
proportion of the algebro-geometric content of the Yang-Baxter equations |33|, ^| . 

One issue for future work is the resolution of the tetrahedron equations, and 
the finding of truly three-dimensional integrability. Indeed there is a competition 
then between the higher overdetermination of the systems of equations to solve, 
and the increasing size of the corresponding group T. This antagonism might be 
resolved either by a relative triviality of the solutions, which are disguised two- 
dimensional solutions, either by a low complexity of the actual realization of T, 
and this leaves room for interesting solutions. The most natural to conjecture is 
that r has a finite realizations, or that its trajectories lie on abelian varieties |79[ . 

The remarkable elliptic parametrization || of the 16-vertex model should 
lead to interesting properties of the corresponding transfer matrices, through 
controlled functional relations on the exact partition function and/or via a 
construction of physically relevant spaces of states. This applies as well to 
three-dimensional models 0, and prompts us to a return to basics, i.e: Bethe 
Ansatz @. 

Acknowledgment. I would like to thank F. Jaeger for a very useful discussion, 
and J. Avan, M. Bellon, J.-M. Maillard, G. Rollet, M. Talon for continued 
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